Rapidly Rotating Bose-Einstein Condensates in and near the Lowest Landau Level 
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We create rapidly rotating Bose-Einstein condensates in the lowest Landau level, by spinning up 
the condensates to rotation rates Q, > 99% of the centrifugal limit for a harmonically trapped gas, 
while reducing the number of atoms. As a consequence, the chemical potential drops below the 
cyclotron energy 2hQ. While in this mean-field quantum Hall regime we still observe an ordered 
vortex lattice, its elastic shear strength is strongly reduced, as evidenced by the observed very low 
frequency of Tkachenko modes. Furthermore, the gas approaches the quasi-two-dimensional limit. 
The associated cross-over from interacting- to ideal-gas behavior along the rotation axis results in 
a shift of the axial breathing mode frequency. 
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Rotating Bose-Einstein condensates (BECs) provide a 
conceptual link between the physics of trapped gases and 
the physics of condensed matter systems such as super- 
fluids, type-II superconductors and quantum HaU effect 
(QHE) materials. In all these systems, striking coun- 
terintuitive effects emerge when an external flux pene- 
trates the sample. For charged particles this flux can be 
provided by a magnetic field, leading to the formation 
of Abrikosov flux line lattices in type-11 superconduc- 
tors [1], or in QHE systems to the formation of corre- 
lated electron-liquids and composite quasiparticles made 
of electrons with attached flux quanta [2] . For neutral su- 
perfluids, the analog to a magnetic field is a rotation of 
the system, which similarly spawns vortices [3] . In rotat- 
ing atomic BECs, the creation of large ordered Abrikosov 
lattices of vortices [4] has recently become possible. 

Here we examine the vortex lattice of harmonically 
trapped BECs approaching the high rotation limit, when 
the centrifugal force quite nearly cancels the radial con- 
fining force. The formal analogy of neutral atoms in this 
limit with electrons in a strong magnetic field has led to 
the prediction that quantum-Hall like properties should 
emerge in rapidly rotating atomic BECs [5] . In particu- 
lar, the single-particle energy states organize into Landau 
levels, and if interactions are weaker than the cyclotron 
energy, primarily the near-degenerate states of the lowest 
Landau level (LLL) are occupied. For rotating bosons 
in the LLL, three regimes have been identified, distin- 
guished by the filling factor i/ = Np/Ny, i.e., the ratio of 
the number of particles [Np) to vortices {Ny). For high 
filling factors (always ^ 500 in our system), the conden- 
sate is in the mean-field quantum Hall regime [6-8] and 
forms an ordered vortex lattice ground state. With de- 
creasing filling factor, the elastic shear strength of the 
vortex lattice decreases, which is reflected in very low 
frequencies of long-wavelength transverse lattice excita- 
tions (Tkachenko oscillations [9-13]). For filling factors 



below V \Q the shear strength is predicted to drop 
sufficiently for quantum fluctuations to melt the vortex 
lattice [5,14], and a variety of strongly correlated vortex 
liquid states similar to those in the Fermionic fractional 
QHE are predicted to appear [5] . For still lower v exotic 
quasiparticle excitations obeying fractional statistics [15] 
are predicted. 

In this Letter we report the observation of rapidly ro- 
tating BECs in the lowest Landau level, and provide ev- 
idence that the elastic shear strength of the vortex lat- 
tice drops substantially as the BEC enters the mean-field 
quantum Hall regime [11]. This effect is a precursor to 
the predicted quantum melting of the lattice at lower 
filling factors. Our rapidly rotating condensates spin 
out into a pancake shape and approach the quasi-two- 
dimensional regime. We observe a corresponding cross- 
over in the spectrum of breathing excitations along the 
axial direction. The high rotation limit has been studied 
experimentally in Ref. [16], focusing on effects of a ro- 
tating trap anisotropy, which is not present in our setup, 
and in Ref. [17] where the addition of a quartic term to 
the trapping potential led to a loss of vortex visibility. 

Our experiments take place in an axially symmetric 
harmonic trap with oscillation frequencies {wp,a;z} = 
27r{8.3, 5.3} Hz. Using an evaporative spin-up technique 
described in Ref. [18] we create condensates containing up 
to 5.5 X 10^ Kb atoms in the |F = l,mp = —1) state, 
rotating about the vertical, z axis, at a rate = 0.95 
[0. = Q,/ijjp is the rotation rate scaled by the cen- 
trifugal rotation limit LOp for a harmonically trapped 
gas). To further approach the limit 17 — )■ 1, wc cm- 
ploy an optical spin-up technique, where the BEC is il- 
luminated uniformly with laser light, and the recoil from 
spontaneously scattered photons removes atoms from the 
condensate. Since the condensate is optically thin to 
the laser light, atoms are removed without position or 
angular-momentum selectivity, such that angular mo- 
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mcntum per particle is unchanged. Atom loss leads to 
a small decrease in cloud radius which, through conser- 
vation of angular momentum, increases O. Over a period 
of up to 2 seconds wc decrease the number of BEC atoms 
by up to a factor of 100, to 5 x 10'*, while increasing Cl 
from 0.95 to more than 0.99 [19]. At this point, further 
reduction in number degrades the quality of images un- 
acceptably. Ongoing evaporation is imposed to retain a 
quasi-pure BEC with no discernible thermal cloud. 

With increasing rotation, centrifugal force distorts the 
cloud into an extremely oblate shape [see Fig.l] and re- 
duces the density significantly - thus the BEC approaches 
the quasi-two-dimensional regime. For the highest rota- 
tion rates we achieve, the chemical potential /x is reduced 
close to the axial oscillator energy, T2D = ~ ^-^^ 
and the gas undergoes a cross-over from interacting- to 
ideal-gas behavior along the axial direction. To probe 
this cross-over, we excite the lowest order axial breath- 
ing mode over a range of rotation rates. For a BEC in the 
axial Thomas-Fermi regime, an axial breathing frequency 
lob = \/3 has been predicted in the limit O — > 1 [20], 
whereas ujb = "i^z is expected for a non-interacting gas. 

To excite the breathing mode, we jump the axial trap 
frequency by 6%, while leaving the radial frequency un- 
changed (within < 0.5%). To extract the axial breathing 
frequency ws, we take 13 nondestructive in-trap images 
of the cloud, perpendicular to the axis of rotation. From 
the oscillation of the axial Thomas- Fermi radius [21] in 
time we obtain ub- Rotation rates are obtained [19] 
from the aspect ratio by averaging over all 13 images to 
eliminate the effect of axial breathing. As shown in Fig. 
2(a), we do indeed observe a frequency cross-over from 
lob = V^ujz to = 2 ^2 as fi — > 1. To quantify under 
which conditions the cross-over occurs, we plot the same 
data vs. [Fig. 2(b)], where the chemical potential is 
determined from the measured atom number, the rota- 
tion rate and the trap frequencies. For T2D < 3, the ra- 
tio lub/^'z starts to deviate from the predicted hydrody- 
namic value, and approaches 2 for our lowest T2D ~ 1-5. 

As ri 1, also the dynamics in the radial plane 
are affected. For the highest rotation rates, interactions 
become sufficiently weak that the chemical potential /U 
drops below the cyclotron energy 2?iO, which is only a 
few percent smaller than the Landau level spacing 2hujp. 
Then, Tlll = 2Hn ^ ^^'^ condensate primarily 
occupies single-particle states in the LLL. These form a 
ladder of near-degenerate states, with a frequency split- 
ting of e = ujp — n. The number of occupied states is 
Nlll ~ We are able to create condensates with 
Flll as low as 0.6, which occupy Nlll ~ 120 states 
with a splitting e < 27r0.06Hz. In this regime of near- 
degenerate single-particle states a drastic decrease of the 
lattice's elastic shear strength takes place. The elastic 
shear modulus C2 is predicted by Baym [11] to decrease 
with increasing rotation rate from its value in the "stiff' 
Thomas-Fermi (TF) limit, Cj-'^ = n(^Q)hfl/8 (where n(o) 



is the BEC number density) to its value in the mean-field 
quantum Hall regime, of C^^^ « 0.16 x Flll x Cj^. We 
directly probe this shear strength by exciting the lowest 
order azimuthally symmetric lattice mode ( (n = 1,to = 
0) Tkachenko mode [9-11]). Its frequency W(i,o) ~ VC2 
is expected to drop by a factor w 2.5 below the TF pre- 
diction when Tlll = 1- 

Our excitation technique for Tkachenko modes has 
been described in Ref. [10]. In brief, we shine a focused 
red detuned laser (A = 850 nm) onto the BEC center, 
along the axis of rotation. This laser draws atoms into 
the center, and Coriolis force diverts the atoms' inward 
motion into the lattice rotation direction. The vortex 
lattice adjusts to this distortion, and after we turn off 
the beam, the lattice elasticity drives oscillations at the 
frequency a'(i.o)- We observe the oscillation by varying 
the wait time after the excitation, and then expanding 
the condensate before imaging the vortex lattice along 
the z-axis [sec Fig. 3(a), (b)]. In Fig. 3(c), we compare 
the measured frequencies W(i,o) to the predictions of Ref. 
[11] for the TF limit and for the mean- field quantum 
Hall regime. For ft < 0.98 {Till > 3), the frequency 
'^(1,0) follows the prediction for the TF regime, whereas 
by O = 0.990 {Tlll = 1-5), W(i g) has dropped to close 
to the prediction for the LLL, thus providing evidence for 
the cross-over to the lower shear modulus C2 predicted 
for the LLL. 

While we are able to produce clouds with Tlll mea- 
sured to be substantially lower than Tlll = 1-5, we are 
unable to accurately measure Tkachenko frequencies un- 
der these extreme conditions, due at least in part to the 
very weakness of C2. The Tkachenko mode frequencies 
become so low that it takes multiple seconds to track 
even a quarter oscillation [Fig. 3(a),(b)], while at the 
same time the very weak shear strength means that even 
minor perturbations to the cloud can cause the lattice to 
melt and the individual cores to lose contrast [22] in a 
matter of seconds. These perturbations can result from 
residual asymmetry of the magnetic trapping potential, 
or from spatial structure in the optical beam used to re- 
duce the atom number, or perhaps from thermal fluctua- 
tions. In contrast, for a "stiff' cloud of 3 x 10^ atoms at 
(l = 0.95 (Tlll ~ 7) we observe that the lattice remains 
ordered, and fl can be kept constant, over the entire 1/e 
lifetime of the BEC (w 3 minutes). 

Finally we discuss observations related to two other 
predicted consequences of rapid rotation. First, the ra- 
dial condensate density profile has been predicted to 
change from the parabolic Thomas-Fermi (TF) profile 
to a Gaussian profile in the LLL [6]. We create clouds 
with varying values of Tlll down to as small as 0.6, 
and after 3 sec equilibration time acquire images of the 
radial density profile. These images fit better to a TF 
profile than to a Gaussian, showing no signs of a cross- 
over in the radial density profile as the LLL is entered. 
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On the theoretical side, it has been pointed out [12] 
that, whereas the predictions of Ref. [6] are based on 
an assumption that the areal density of vortices is uni- 
form, a very small deviation from uniformity could ac- 
count for an overall TF density profile. Furthermore, 
according to Ref. [8], a Gaussian radial profile is to be 
expected only if Na/R^ <C 1 {N: number of atoms in 
the BEG, a: scattering length, R^: axial BEG radius), 
whereas our measurements were performed on conden- 
sates with Na/Rz ^ 26. It certainly seems reason- 
able to expect a TF profile for a given direction when 
the chemical potential is much greater than the confin- 
ing frequency in that direction, and for our condensates 
with Tlll = 0.6, the chemical potential fi/H « 27rlOHz 
while the (centrifugally weakened) effective radial trap 

frequency We// = y'w^^-Q^ Ri 27r 1 Hz. 

Second, there has been a debate about the possibility 
of vortex core overlap in the high rotation limit [7,8,23]. 
To achieve adequate imaging resolution, we expand the 
BEG before imaging the vortex lattice. Our expansion 
procedure [18] leads to nearly pure two-dimensional (2- 
d) expansion (a thirteen-fold expansion in radius occurs 
while the axial size changes by < 25%). We expect that 
under these conditions, the ratio of the vortex core area 
to the lattice's unit cell area should be preserved through- 
out the expansion process. We fit a 2-d Gaussian to the 
missing optical density associated with each vortex core. 
The vortex radius is defined to be the RMS radius 
of the 2-d Gaussian, which is 0.60 times its FWHM. We 
then define the fractional area A occupied by the vor- 
tices to be ^ = nyirrl, where n„ is the areal density of 
vortices. In the limit of many vortices, the theoretical 
prediction for n„ is mfl/{iTh). 

To determine a theoretical value for the vortex core 
size, we perform a mimerical simulation of the Gross- 
Pitacvskii equation of a BEG containing an isolated vor- 
tex. Wc obtain r„ = 1.94 x ^ with ^ = {8nna)^2^ 
where a is the scattering length and n is the density 
[24]. The resulting prediction for A can be expressed 
as A = 1.34 X {Tlll}^^- This value exceeds unity for 
Tlll < 1.34, which has led to the prediction that vor- 
tices should merge in the LLL limit. An alternate treat- 
ment due to Baym and coworkers [8] on the other hand 
predicts that A saturates at 0.225 in the LLL. Our data 
for A are plotted in Fig. 4. For F^^^ < 0.1 the data agree 
reasonably well with our numerical result. For larger 
^LLL data clearly show saturation of ^ at a value 
close to the LLL limit [8] , rather than a divergence of A 
as fi ^ 1. Further details on vortex core structure will 
be provided in a future publication [25]. 

In conclusion, we have created rapidly rotating BEGs 
in the lowest Landau level. The vortex lattice remains 
ordered, but its elastic shear strength is drastically re- 
duced. In expansion- images we find no divergence of 
vortex core area as f2 — > 1, as well as no deviation from 



a radial Thomas-Fermi profile. Additionally, our rapidly 
rotating BEGs approach the quasi-two-dimensional limit. 
We remain far from the regime in which quantum fluc- 
tuations [14] should destroy the lattice, but observing 
the effects of thermal fluctuations [26,27] in this reduced- 
dimensionality system may be possible. 
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FIG. 1. Side view images of BECs in trap, (a) Static BEC. 

The aspect ratio Rz/Rp = 1.57 (A'' = 3.8 x 10*^ atoms) resem- 
bles the prolate trap shape, (b) After evaporative spin-up, 
N = 3.3 X 10®, Q, = 0.953, (c) evaporative plus optical spin-up, 
N = 1.9 xlO^,n = 0.993. Due to centrifugal distortion the 
aspect ratio is changed by a factor 8 compared to (a). 
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FIG. 2. Measured axial breathing frequency wb/cj^ (a) as 
a function of rotation rate f2 and (b) vs. T2d- Solid line: 
Prediction for the hydrodynamic regime [20]; Dashed line: 
ideal gas limit. For tl > 0.98 (r2D < 3) a cross-over from 
interacting- to ideal-gas behavior is observed. Representative 
error bars are shown for two data points. 
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FIG. 3. (a). 



(b) Tkachcnko mode at Tlll = 1.2 
{N = 1.5 X lO'', Cl = 0.989): (a) directly after excitation, 
(b) after 1 sec - the lattice oscillation has not yet completed 
1/4 cycle, (c) Comparison of measured Tkachenko mode fre- 
quency ct;(i.o) (solid symbols) vs. Q, to theory [11], using vortex 
lattice shear modulus Cj^ in the Thomas-Fermi (TF) limit 
(circles), and C|^^^ in the mean- field quantum Hall regime 
(stars). Note that both N and Tlll decrease as Q increases. 
For Flll ~ 3 (reached at = 7.8 x 10\ O. ^ 0.978) the data 
cross over from the TF to the quantum Hall prediction. 
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FIG. 4. Fraction of the condensate surface area occupied 
by vortex cores, A, measured after condensate expansion 
(squares), plotted vs. the inverse of the lowest Landau level 
parameter, {Tlll)~^ = 2hfl/n. The data clearly show a sat- 
uration of ^, as — > 1. Dashed line: prediction (see text) 
for the pre-expansion value at low rotation rate. Dotted line: 
result of Ref.[8] for the saturated value of A in the LLL. 
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